Wigner Crystallization in inhomogeneous one dimensional wires 
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We explore the theory of electrons confined by one dimensional power law potentials. We cal- 
culate the density profile in the high density electron gas, the low density Wigner crystal, and the 
intermediate regime. We extract the momentum space wavefunction of the electron at the Fermi 
surface, which can be measured in experiments on tunneling between parallel wires. The onset of 
localization leads to a dramatic broadening of the momentum space wavefunction together with 
pronounced sharpening (in energy) of the tunneling spectrum. 
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Advances in microprocessing have enabled scientists to 
construct ultra-high mobility one dimensional wires with 
transverse size d±_ ~ 20nm, which is sufficiently small to 
freeze out all transverse electronic motioni . By allowing 
tunneling between two such parallel wires, Auslaender et 
al^ have been able to map out the spectrum of elemen- 
tary excitations, providing good evidence of spin-charge 
separation. Further experiments have found evidence of 
electronic localization in these wires when a gate elec- 
trode is used to deplete the density of electrons^. This 
localization is possibly caused by an interaction driven 
transition into a Wigner crystal stated, where each elec- 
tron is localized by repulsion from its neighbors. Here 
we investigate the theory of such crystallization in one 
dimensional wires. 

The wires used in these experiments are not infinitely 
long. Their finite length, which is controlled by electro- 
static potentials, leads to striking fringes in the tunneling 
spectrumi. These fringes cannot be explained by parti- 
cles bounded by 'infinite walls', but require a more ac- 
curate modeling of the confinement. Tserkovnyak et ali 
found that the experimental results were consistent with 
a power-law potential V ex t(x) = Vox 13 , with 6 < (i < 7. 
In discussing experimental signatures of crystallization 
in these systems, we must take into account this same 
physics, and explore the interplay between the external 
confining potential and the interparticle Coulomb repul- 
sion. Due to the influence of the gate electrode, the elec- 
trons in the localization experiment*^ feel a complicated 
potential, with two minima separated by a barrier. For 
most of our analysis we use the simpler power law poten- 
tial. 

We are interested in gaining an intuitive understanding 
the crossover between electron liquids and solids in these 
inhomogeneous wires. We therefore choose to work with 
the simplest techniques which can capture the qualitative 
features of these experiments. We use a Thomas-Fermi 
description of the high density electron gas (sec. IHIfl : a 
semiclassical description of the Wigner crystal (sec. Ej i 
and a Hartree-Fock description of the cross-over (sec. Iv)l . 
These models provide direct insight into the experimental 
signatures of this transition. They do not, however, give 
direct information about excitations, transport proper- 
ties, or fractionalization of the electrons into spin and 



charge degrees of freedom. 

Our approach complements prior theoretical work, 
which predominantly treated homogeneous infinite 
systems^, small numbers of electrons^, or infinite square 
well potentials^*. Particular note should be taken of the 
recent work by Fiete et al^i which analyzes the experi- 
ments discussed here*** in the limits of extremely small 
and extremely large particle numbers. 

Our main accomplishment is that we are able to fol- 
low the evolution of the quantum state from the weakly 
interacting high density regime to the strongly interact- 
ing low density regime. We describe the electronic den- 
sity profile, and see how Friedel oscillations evolve into 
Wigner-crystal correlations. In all of these regimes, we 
predict the tunnelling spectrum between parallel wires, 
giving qualitative explanation of the striking features in 
recent experiments^**. 



I. MODEL 

We consider a gas of electrons with a = 2 spin 
states confined to a one dimensional wire, experiencing 
an external potential V ox t(x), where x is the coordinate 
along the wire. The electrons interact through pairwise 
three dimensional coulomb interactions, Uy,d{fi ~ r j) — 
e 2 /(47re|ri — r.j\) — (h 2 /m*a)l/\ri — rj\, where e is the 
dielectric constant of the medium and m* is the elec- 
trons effective mass. Using parameters for n-doped GaAs 
(e ~ 13, m* w 0.067m e , where m e is the free-space 
electron mass), the effective Bohr radius is quite large, 
a w 10~ 8 m. In this wire, the electronic wavefunction 
has a transverse size d± ~ 20nm. Integrating out these 
transverse dimensions, the electrons feel a regularized 
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Coulomb interaction U(r) = (h /m*a)f ( i ± (ri 

useful regularizations are f^\x) = (x 2 + c 

fj^\x) — Mindxl^ 1 , d -1 ), where Min(x, y) is the smaller 
of x and y. As d± — > the exact form of the regular- 
ization becomes unimportant. Here we mostly rely on 

One estimates the importance of correlations caused 
by the Coulomb interactions by comparing the interac- 
tion energy between two neighboring particles E- m t = 



2 



h 2 n/am* to the Fermi (kinetic) energy £/ = k 2 /2m* — 

2TT 2 H 2 n 2 /am* , where n is the one dimensional density 
of particles, and we have assumed nd± <C 1. Kinetic 
energy dominates at high densities, when na ^> l/2ir 2 . 
In that limit one expects to find a liquid state, where 
electrons are delocalized. Conversely, when na <C l/2n 2 , 
the Coulomb interaction dominates, and a Wigner crys- 
tal should be formed. 

As a compromise between experimental relevance 
and simplicity we concentrate on power law potentials 
Vext(x) = (h 2 /2m*)xPw~P~ 2 , where w, which parame- 
terizes the potential strength, roughly coincides with the 
size of the single particle ground state in this potential. 
Figure ^b) illustrates the structure of this potential for 
different values of (3. More complicated potentials are 
briefly discussed in appendix 151 

For the relatively flat (large (3) potentials used in ex- 
periments, w also roughly coincides with the 'lithographic 
length,' Liith, which is the physical distance between the 
gates which generate the potential barriers. The limit 
[3 — > oo corresponds to an infinite square well. Anal- 
ysis of experimental data by Tserkovnyak et ali show 
that w ~ 3/im. Taking the density of electrons to be 
n w N/w, one expects crystallization when N < N*, 
with N* ~ 20. 



II. TUNNELING BETWEEN PARALLEL WIRES 

The experiments of Auslaender et al^ measure the tun- 
neling current between the short wire described in sec- 
tion [I] and a parallel long wire, which we will take to 
be infinite and uniform. Momentum is conserved in the 
tunneling. 

We will make the following approximations to analyze 
this experiment. (1) We will treat the tunneling within 
perturbation theory. This is an important limit to study, 
though, as will be shown below, the experiments are in 
a regime where higher order processes may be playing a 
role. (2) We will neglect correlations between electrons in 
the two wires. One expects that this approximation will 
break down in the extremely low density limt, when the 
separation between electrons in one of the wires becomes 
comparable to the interwire spacing. (3) In the bulk 
of this paper we use an independent electron picture to 
study the electrons in the short wire. The consequences 
and validity of this approximation are discussed in sec- 
tion ^TC] and are explored in more detail in appendix lAl 
We expect that in the limit where the splitting between 
different spin states is small compared to the tempera- 
ture or linewidths that the independent electron picture 
correctly predicts the magnetic field dependence of the 
tunneling current. However, this approximation results 
in a vast overestimation of the overall amplitude of the 
current. 

We consider the geometry sketched in figure ^ The 
wires are separated by dt ~ 6nm, and the barrier sep- 
arating the wires has height U ~ 300meV, arising from 



the conduction band offset between the GaAs wires, and 
the intervening AlGaAs. Using gate electrodes, the ex- 
perimentalists control the chemical potential difference 
SV between the wires. 12 A magnetic field, perpendicular 
to the wires, gives a momentum kick Q = eBd t /Ti to a 
tunneling electron^. We therefore consider the tunneling 
Hamiltonian, 

Ht^-T^J ^ sv m(k + Q)tpa(k) + H.C., (1) 

where Vv(fc) — f dx e~ lkx ip a (x) and (j> a (k) 
j dxe~ tkx (f> a (x) are respectively the operators which an- 
nihilate particles with momentum k and spin a in the 
short and long wire. The magnitude of the tunneling 
matrix element T is estimated by examining the energy 
states in the double square well geometry sketched in fig- 
ure n This one-dimensional potential is a crude model 
of the physics transverse to the wires. We can identify 
T « SE/2, where 8E is the energy splitting between the 
two lowest energy single particle states. In the limit of 
a deep, wide barrier (U > (h 2 / 'm* d 2 ) , (H 2 / 'm* d\)) , one 
finds 

6E/2 = ^-e-^, (2) 
m*d ± K 

where h 2 n 2 /2m* — U (so numerically k w 0.2/nm). 
which gives T « 30meV, which should be compared 
to the spacing of states transverse to the wire, SE = 
4ir 2 h 2 /m*d 2 L w fOOmeV. Due to this separation of scales, 
we treat H t perturbatively. Note, however, that the ra- 
tio T/SE is not extremely small, so that quantitative 
comparison with experiment would require going beyond 
lowest order in perturbation theory. 
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FIG. 1: (Color Online) Model of potential, (a) Two dimen- 
sional grayscale plot, where darker colors represent lower po- 
tential energy, (b) Potential along one wire [dashed blue line 
in (a)] for j3 = 2, 3, 6, 10, 50, 1000. (c) Potential transverse to 
the two wires [dotted red line in (b)]. 
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A. Formal expression for tunneling current 

The current operator, denned by I(t) = —dN^/dt, 
where = J" V^(fc)Vv(fc) is the number of par- 

ticles in the short wire, is given by 

dk 



I(t) = ~T / ^E Im e uvt tUk + Q,t)Mk,t) 



To lowest order in T, the tunneling current is 



(I{t)) 



1 



dr 7 e(t-T){[I(t),H t (T)}) 



T 2 Ai(Q,uj = 5V), 



(3) 



(4) 
(5) 



tion by Aj(q,uj) = 2Imxf(q,oj), where = 
/die*"(0(t)/i)[xf((i,t) - xf(MM with X >{<l,t) = 
/# ^EM(k + q,t)Mk,t)ti(k>,0)Mk> + q,0)) and 
Xf(q,t) = J^^EM(k',0)Mk' + q,0)$t(k + 
q,t)ijja-{k,i)}. To simplify this expression we assume that 
in the electrons in the two wires are independent. Thus 
we neglect correlations between the two wires caused by 
the Coulomb interaction. Such correlations are likely to 
be present in the extreme Wigner crystal limit, where the 
spacing between electrons in one of the wires is large com- 
pared to the wire separation. Quantitative comparison 
with experiments in this limit would require including 
these correlations, which will partially screen the long- 
range part of the Coulomb interaction. 



where the spectral density of current fluctuations 
is related to the retarded current response func- 



Within our approximation the expectation values can 
be factored, leading to 



Ai{q,w) 



dk dk 1 
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G< a {k, k'- y)G> a {k + q ,W + q;v-u)) -G> a (k, k'; v)G<„{k + q, k' + q; v - u>) (6) 



[/(f) - f( v - w )] Ap<r(k, k'; v)A^ a {k + q, k' + q; v - u), 



(7) 



where the single particle Greens functions 
G>/< (*:,*» = jdte^G^ik.k'-t), are given 
by G> CT (fc,fc';t) = (^(A,t)4(fe',0)), = 
— (<j)v(k',0)(f) a (k,t)}, and equivalent expressions hold 

for G^^(k,k' ,uj). These Greens functions are related 
to the appropriate single particle spectral density by 
G>(lu) = [1 - f(u)]A(u) and G<(u) = f(u)A(u), where 
f{ijj) = [e 13 ^-^ + l}- 1 is the Fermi function. 

If we use the free electron spectral density for the long 
wire, A(p a (/c, k'\ w) = (27r) 2 c5(fc - k')5{u) - k 2 /2m*), one 
arrives at the simple result that the tunneling current is 
a direct measure of the single particle spectral density in 
the short wire, 

AMq) = |^^[/(c + fc 2 /2m*)-/(fc 2 /2m*)] 

xA^ )(r (k-q,uj + k 2 /2m*) (8) 

/dk 
— Y^fk^-fk^^k-q^/Zm*), 



r 



where the last line neglects terms of order uj 2 , and uses 
f k = f(k 2 /2m*). 

At zero temperature, the Fermi functions become step 
functions and to lowest order in u, equation (jSJ becomes 



V f rr 



+A i)(T (-^2m*£f-q,Sf) 



(9) 



B. General Features 



From the definition A^, a -(k,uj) 



G5 cr (fe, cj), the single particle spectral density can be 
written as 



A^ a {k, lo) 



2tt 
~Z 



if 



P^Ei-^Ni) 



(i\Mk)\f) 8{Lo-{E f -E t ))~ (f\Mk)\i) S(uj - {E % - E f )) 



(10) 



where \i) represents a many body state containing 
Ni electrons and possessing energy Ei. Normaliza- 



tion is given by the grand partition function Z 



4 



•^2^ e ~(i(Ei-^Ni) ^ zer0 temperature, the sum over i 
is omitted, and \i) is replace by the ground state. 

One thus sees that current flows whenever momentum 
and energy can be conserved in the tunneling process. 
The magnitude of the current is set by the degree of 
overlap between states with different particle number. 

When the electrons become localized in a Wigner crys- 
tal, the overlap between states with different particle 
number becomes extremely small (because removing a 
particle results in rearrangement of the entire crystal - 
a form of orthogonality catastrophe). Consequently, the 
average tunneling current should drop. This drop in aver- 
age current is accompanied by a sharpening of line- widths 
and does not lead to a complete loss of the experimen- 
tal signal. Instead, it leads to "coulomb-blockade" type 
peaks, where current only flows for discrete values of SV . 

One will see discrete peaks if the line width is small 
compared to the level spacing in the short wire, 6E ~ 
(2ir) 2 h 2 /m*w 2 ~ 0.05meV. Assuming that the lifetime is 
only due to coupling between the two wires, then the line 
width is given by an expression very similar to that for 
the tunneling current. Since T is large on the scale of 5E 1 
one will only see these discrete peaks once the overlaps 
in I|1U|I become extremely small. 



C. Independent particle approximation 

In an independent particle picture (such as Hartree- 
Fock), A^ a (k,u>) = £\ *{k)\ 2 2tt5{uj - E j<7 ), where 
4>j,<r{k) is the momentum space wavefunction of the j's 
single electron orbital of spin a particles, with energy 
Ej (j . In general <j)j j ^ 4>j i, and the orbitals are cho- 
sen self-consistently. Within such an approximation the 
tunneling current is 

I = I J2\^a{k)\ 2 6(E na -£f), (11) 

I la 

where, I = T 2 SV/{2n) v / m*/2£ f , k = y/2m*£ f -Q, and 
we have assumed that Q = eBdt/h > 0. The tunneling 
current therefore measures the momentum space density 
of the particles at the Fermi surface. If there are no 
particles at the Fermi surface, then no current flows. 

The delta functions in (|llfl are broadened by finite 
lifetime and finite temperature. In particular if there 
are two different states i,j with energies Ei/j, and in- 
verse lifetimes T, which satisfy \Ei/j — £f\ <C T or 
\E i /j — £f \ <C ksT, then the current is proportional to 
the incoherent sum 

/ex |^(fc)| 2 + |^(fc)| 2 . (12) 

This feature will be important in high symmetry situa- 
tions where there is a near degeneracy (see, for example, 
section IV Af) . 

In the remaining sections of this paper, we calculate 
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FIG. 2: Interaction induced change of dimensionless density 
n of spin- 1/2 electron gas confined to a one-dimensional po- 
tential V(x) oc x 6 . Interaction parameters are a = 10 and 
d — 10~ 3 . Solid line: Equation 115H using regularization /' 2 '; 
dashed line: approximation described above eq. H16I : dotted 
line: noninteracting profile. Inset: Renormalized dimension- 
less cloud radius i as a function of 1/a from solving eq. HIGH . 

In its conventional form, the independent particle ap- 
proximation does not include the effects of orthogonality 
between many-body states of different particle number, 
as described in section Til Bl In appendix 1X1 we explore a 
simple example which lets us systematically study these 
effects. We find that the omitted physics predominately 
modifies the overall scale of the tunneling current, and 
can be partially accounted for by renormalizing Iq in 
equation l|ll|l . Hartree-Fock also misses fine-detail fringe 
structure in the magnetic field dependance of the tun- 
neling current. The prediction of the envelope of the 
tunneling current appears to be qualitatively correct. 



III. HIGH DENSITY: THOMAS-FERMI 

The simplest picture of a gas of electrons in a confined 
geometry comes from the Thomas-Fermi (local density) 
approximation, where the system is described by a lo- 
cal chemical potential ^jl{x) = /xo — V c s{x). The global 
chemical potential is no = fi 2 k 2 /2m*, and the effective 
potential V c s(x) includes both the external field and the 
interactions between the particles. Within the Hartree 
approximation, which is valid at high densities no > 1, 
the effective potential is 

V c ff(x) = J ' dx ' U(x — x)n(x). (13) 

The chemical potential is then related to the density 
by the relationship for a homogeneous gas, /i(r) = 
2ir 2 h 2 n 2 /a 2 m* . Self-consistency requires that the den- 
sity obeys a nonlinear integral equation, 

W) = U»)+ fdx'U(x-x')n(r') + ^^n 2 (r). (14) 
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This local density approximation only makes sense if 
the density changes slowly compared to the interpar- 
ticle spacing [(d x n(x)) / n{x) 2 <C 1]. We write lTH|l 
in dimensionless form by introducing k 2 — 2m* /iq/H 2 , 
R 13 = k 2 w 2+l3 , f = r/R, a = k a a, n = k Q n, d = d/R, 

2 



47T" 



n 2 {r) = 1 - 



f d {f-f')n{f')df'. (15) 



At fixed chemical potential the interactions can only 
reduce the density. Therefore the density is al- 
ways bounded by the noninteracting result, uq — 
(cr/27r)vl — f@ . In particular, the density always van- 
ishes for f > 1. For very steep potentials, (3 — > oo, the 
noninteracting Thomas-Fermi radius R approaches the 
trap length, R/w « 1 + 2/T 1 log(fc «;) + O^ 2 ). 

For a ^> 1, equation 115|) can be solved iterativeljiii 
In this limit, the interactions predominately renormalizc 
the chemical potential. Approximating the integral in 
(|15|l by its value at f = 0, one self-consistently finds that 
(2i:/a)n(f) w \Jt 2 — W where t solves the transcendental 
equation 



i 



: t(A + 41ogt). 



(16) 



The parameter A depends logarithmically on the cut- 
off d±. For the regularization f%, A = 1 — 4//3 + 
(2//3) log(4/dj_). Figure [3 compares this approximation 
to the numerical solution of l|15[l . and the inset shows t 
as a function of a. 



A. Wavefunction of last occupied state 

As previously explained, within an independent par- 
ticle picture the tunneling current between two parallel 
wires is proportional to the momentum space wavefunc- 
tion of the highest occupied single particle state. This 
momentum space wavefunction will be strongly peaked 
about kf(r = 0), which can be approximated as k^t, 
where h 2 k 2 /2m* = fi , and t is given by eq. JTHJ. Fig- 
ure shows kf = kf (r = 0) as a function of [Iq for rea- 
sonable parameters. 

To go beyond this single wavevector approximation, 
we use semiclassical means to calculate the wavefunction 
of the state at the Fermi surface. We find the modulus 
by imagining that we reduce the chemical potential by a 
very small amount, so that this last state goes from filled 
to empty. The density change must coincide with the 
density n,f(r) of that last single particle state, implying 
that 



n f (x) = \i> f (x)\ 2 = n 



dn(r) 
d[i 



(17) 



where is a normalization constant determined by set- 
ting Jnf(x)dx = 1. This wavefunction should have a 
phase whose derivative gives the local fermi momentum. 



dx 




FIG. 3: Peak wavevector kf — tko of wavefunction at fermi 
surface as a function of chemical potential = ft 2 fco/2m* 
within the Thomas-Fermi approximation. Dotted line shows 
noninteracting result, t = 1. Solid line includes Coulomb 
interactions with a — lOnm, d± = 20nm, in a power law trap 
with /3 = 6 and w = 3/xm. This approximation will break 
down when kf ^ 1/a. 



The + and — solutions correspond to left-moving and 
right moving waves. In principle, one must combine 
these two solutions to form a standing wave. Within 
the Hartree approximation, this procedure is equivalent 
to finding the WKB wavefunction in the self-consistent 
potential. We do not discuss the details as when we com- 
pute the momentum space wavefunction, only the struc- 
ture near k = will be affected by how we superimpose 
these two solutions. 

In the case where the Coulomb interaction can be ne- 
glected, equation i|17fl reduces to \i[)f(x)\ cx kf(x)~ 1 ^ 2 cx 
Wmo — V ex t(x), and equation lfT%|l becomes dip/dx — 
kf(x). Not surprisingly, in this limit ipf(x) is exactly 
the WKB wavefunction for an electron with energy /j,q 
in the potential V ex ±{x). Our arguments therefore re- 
duce to the semiclassical arguments of Tserkovnyak et 
ali for the noninteracting gas. The analysis is very 
similar for our approximation where n 2 cx t 2 — f 3 . 
Since interactions only renormalize the chemical poten- 
tial, the derivative dn 2 (r)/d^Lo is a constant, indepen- 
dent of space, and we can use the noninteracting result 
with a renormalized chemical potential. That is, we take 
k f (x) = yjt 2 nv - V cxt {x) = k Vt 2 - xP. 

Following Tserkovnyak et al»Z>, we calculate the mo- 
mentum space wavefunction, <j)f{q) = fdxipf(x)e~ lqx by 
Laplace's method. The Fourier integral is dominated 
by regions of space near where q = ±fcj(x). Moreover, 
since <fif(q) is peaked around q = ±fc/(0), we can ex- 
pand <p(x) near x = 0. Since |</>/(<?)| is symmetric, there 
is no loss of generality in taking q > 0. We introduce 
p = {kf(Q)-q)X, with XP+ 1 = 2k f w l3+2 . Changing vari- 
ables to y = px/X, and assuming (fc/(0) — q) <C fe/(0), 
we find 



±-\/2m* n(x) = kf(x) 



(18) 



9.' 



= Re 



dyexp ipy - i-r-— 
o V P + 1 



(19) 
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where il' is a normalization constant. In the case (3 — 2, 
the integral in l|19fl is an Airy function, though for other 
values of it is not a familiar special function. The 
argument of the exponential, x, is stationary at y = yo = 
pi/0. The curvature is x"{ya) = — i0p 1 ^ 1 ^ 13 - This is an 
isolated saddle point when y 2 \x"(yo)\ — 0\p\ 1+1 ^ ~> 1, 
in which case 



4>f{q) 



Re 



2tt 



i0p 1 - 1 /f } 



exp 



1 + 



(20) 



By deforming the contour of integration, this equation 
holds for both p > and p < (in the latter case one 
should take the principle branch, giving p a small nega- 
tive imaginary part). From l|20l) . one sees that if ^ oo, 
the wave function falls off exponentially as p — > — oo. 
For p <C fiP/( 1+l3 \ instead of expanding about the saddle 
point, we expand about y — 0, finding as p — > 0, 



0' 



= Re 



(/3 + 1) -/V(/m) r (_l_ 



(21) 



e 



-t(7T/2)(l/(/3+l)). 



IV. LOW DENSITY: WIGNER CRYSTAL 

In section ITITl we neglected the fact that electrons are 
discrete entities and that an electron does not interact 
with itself. Including this discreteness, the interaction 
energy can be reduced by localizing each electron to a 
small volume - maximizing the distance between the 
charge distribution of neighboring particles. Since such 
localization costs kinetic energy, it is only favorable in 
the low density gas where interaction energy dominates 
over kinetic energy 

To discuss this phase, we imagine that we have N elec- 
trons, localized at positions Xj where j = 1,2, ... ,N. Ne- 
glecting kinetic terms, the energy of a given configuration 
is 



(1/2) $>( 



(22) 



As with the profile in the dense limit (fig. |5J), we nu- 
merically find that the resulting density is quite flat, 
and is well approximated by evenly spaced electrons, 
Xjn w j - (N + 1) /2, with uniform density n. Minimizing 
(1221) with respect to n yields 



2w 

IT 

2w 
~N 



+ l w 
20 ~a 
1 



N 2 [ logN- i 



i/OM-i) 



(23) 
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FIG. 4: Interaction dependence of the momentum space wave- 
function of last single-particle state within the Hartree-Fock 
approximation. Eighteen particles are confined to the wire. 
Lighter colors correspond to larger values of |<^(fc)| 2 . The 
horizontal axis shows different interaction strengths parame- 
terized by the ratio of the Bohr radius a to the wire length 
w (note the logarithmic scale). The potential has exponent 
j3 = 6, and a cutoff, d — 0.07™. 



In this Wigner crystal state, each electron is confined 
to a length £ < n^ 1 . To calculate this length, and the 
wavefunction of the localized electron, we expand the po- 
tential in eq. (1221) for small fluctuations in the position 
of one of the particles, 



dx 2 oxt 



2h 2 



ma 



(24) 



We approximate the sum over a finite number of particles 
by that of an infinite chain, taking \%i — %j\ ~ 

2n 3 J2'iZi I 3 — 2£(3)n 3 . Numerically, the Riemann zeta 
function is £(3) » 1.2. The curvature of the external po- 
tential is generically negligible in this limit, hence each 
electron is trapped by a harmonic oscillator potential 
with m*w^, c = (h 2 j m* a)4£(3)n 3 and has a wavefunction 
ipj(x) oc exp(-(x - Xj) 2 /2£ 2 ) with £ 4 = (h/m*uj wc ) 2 = 
(m~ 3 /(4C(3)). This implies that the momentum space 
wavefunction has modulus \tpj(k)\ 2 oc exp(— k 2 £ 2 ). 

Comparing this momentum space wavefunction to the 
high density predictions [equation (|19fl through JUJ], we 
see that in the Wigner crystal, | </>(<?) | 2 1S more spread out 
in momentum space, and does not contain any oscilla- 
tions. 



V. HARTREE-FOCK STUDY OF CROSSOVER 

We investigate the crossover between the high density 
electron gas and the low density Wigner crystal within 
the Hartree-Fock approximation. In the extreme lim- 
its, Hartree-Fock reduces to our previous approxima- 
tions. It also provides semi-quantitative understanding 
of the intermediate regime. In uniform higher dimen- 
sional systems, Hartree-Fock overestimates the stability 
of the Wigner crystal^, and one therefore expects some 
systematic errors as compared to an exact many-body 
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FIG. 5: (Color Online) Real space electron density n. Param- 
eters coincide with Fig.0] Solid, blue dotted, and red dashed 
lines correspond to total density and the density of up and 
down spin electrons [In (a), the latter two curves coincide]. 
Panels (a), (b), (c) correspond to a/w— 0.5, 0.15, 0.05. 



calculation. Unlike bosonization and Luttinger Liquid 
techniques, Hartree-Fock correctly accounts for band cur- 
vature. Band curvature plays a significant role at low 
densities. Appendix 1X1 compares Hartree-Fock and exact 
results for two particles. 

Hartree-Fock is most simply thought of as a varia- 
tional method, where one searches for the Slater deter- 
minant which minimizes the energyi^. This procedure 
is highly nontrivial as the energy landscape in the space 
of Slater determinants is quite complex with many local 
minima. We start by considering the limit of vanishing 
interactions (a/w — > oo), where we find the exact solu- 
tion by discretizing space and solving the single-particle 
Schrodinger equation. With fixed numbers of up and 
down spin electrons (rif and n.jj, we then gradually de- 
crease w / a. For each value oiw/a we iteratively solve the 
Hartree-Fock equations in discretized space. Although 
there is no guarantee that the states we find in this man- 
ner are absolute energy minima, they consistently have 
lower energies than all other states that we have found 
by iterating the Hartree-Fock equations from different 
starting points. 

An advantage of the Hartree-Fock approximation is 
that since it involves an independent electron approxi- 
mation, one has direct access to single-particle observ- 
ables, such as the wavefunction of the last bound state. 
Figure0]shows the momentum-space wavefunction of the 
last bound state for a system of eighteen particles (nine in 
each spin state). For large a/w, this last single-particle 
state is delocalized in real space, resulting in a series 
of momentum-space peaks which are well described by 
equation (|19fl . For small a/w, this last state is localized, 
resulting in a broad spread of momenta. Illustrative real- 



space density profiles are shown in figures |S] a, b,c. In (a) 
we see that even for arbitrarily weak interactions the real 
space density is corrugated. The corrugations are caused 
by the free electron response to an inhomogeneous poten- 
tial and are analogous to the Friedeii& oscillations seen 
in electron density near an impurity. The asymptotic be- 
havior of these oscillations has been analytically calcu- 
lated for harmonic confinemenljii. There are nine peaks 
corresponding to the nine doubly occupied single particle 
states. In (b) we see that as we increase the interaction 
strength each of the density peaks split into two, so there 
is one peak per particle. A spin-density wave appears. 
In (c) we see a Wigner crystal with an antifcrromagnetic 
spin profile. The density varies smoothly as one increases 
the interactions. 



A. Fringes 

When an odd number of particles form a Wigner crys- 
tal the highest occupied single particle state, with wave- 
function <pi(r), is in a symmetric superposition of two 
locations, as illustrated in figure and has fringes in its 
momentum space wavefunction, cf>i(k). The next highest 
energy state, 4>j(r), is in the antisymmetric superposition 
of these same two locations, and has momentum space 
fringes which are 90 degrees out of phase with those for 
4>i(k). Since the splitting between these two states is ex- 
ponentially small, the tunneling current is proportional 
to |0i(fc)| 2 + |0j(fc)| 2 (see section Hi C|) . which, as illus- 
trated in figure|nt, contains no fringes. 

Going beyond Hartree-Fock, one generically expects 
the many-body wavefunction to include amplitudes for 
each particle to be in a superposition of more than one 
place. Thus, in an exact treatment one would always ex- 
pect that tunneling from any one particular state to any 
other particular state would result in fringes. However, 
as with the Hartree-Fock result, if one incoherently adds 
the fringes from several symmetry related states then the 
fringes will be washed out. In particular, if the tempera- 
ture or lifetime is greater than the splitting between these 
states then no fringes will be observable. 



VI. SUMMARY 

We have presented simple analytic theories of very high 
and very low density electrons in an inhomogeneous one 
dimensional potential. At high density an electron liquid 
is formed where the density displays Friedel oscillations 
from response to the potential [Fig. Ufa)]. At low den- 
sities a crystal-like state is formed where the electrons 
are well-separated [Fig. Eta)]. Using the Hartree-Fock 
approximation we numerically studied the cross-over be- 
tween these limits. Intermediate between the two ex- 
tremes we found a spin-density- wave state [Fig. E£b)] . 

In the extreme high density limit our approach is exact. 
In the extreme low density limit, our approach will give 
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Although simple mean-field analysis (sec. IV A|) and 
more general symmetry arguments^ suggest that the 
large peak at k = will show a fringe structure, we 
have demonstrated that if either the temperature or the 
inverse state lifetime T is larger that the splitting be- 
tween spin states then this structure is washed out, and 
may therefore be extremely hard to observe experimen- 
tally. We should emphasize that the Hartree-Fock theory 
used here is unable to accurately describe these fringes. 
From our studies of the two-particle system, we believe 
that Hartree-Fock does describe the envelope of the peak, 
though this needs to be confirmed by studying larger sys- 
tems. 
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FIG. 6: (Color Online) (a) Density of 13 electrons with a/w = 
0.1, see caption of fig|S]for key. (b) Wavefunction <f>i(x) (solid) 
and <f)j( x ) (dashed) of the two highest highest energy occu- 
pied states, (c) Momentum space wavefunctions: red dashed, 
\Mk)\ 2 ; blue dotted, |<fe(fc)| 2 ; solid, (j<Mfc)j 2 + \4>j{k)\ 2 ). 



the exact electronic density. Our mean-field approach, 
however, can not describe the quantum spin correlations 
which occur in the Wigner crystal state. Instead, it pre- 
dicts a classical Neel order. One expects that our descrip- 
tion of the crossover is qualitatively correct, but that the 
stability of the crystalline state is overestimated. 

Within the independent electron picture we have calcu- 
lated the tunneling current between electrons in a short 
and long wire. In the liquid phase, the zero bias momen- 
tum resolved tunneling is dominated by sharp peaks at 
±fc/, and there is large overlap between states with differ- 
ent particle number. As interactions increase, these two 
peaks broaden and merge. In the crystaline phase, the 
tunneling displays a large feature centered at k = 0. This 
feature is attributed to electron localization. Going be- 
yond our mean-field picture, it is clear that in this limit 
the overlap between states with different particle num- 
ber becomes extremely small, which leads to coulomb- 
blockade structures. 

For these tunneling predictions we work in linear re- 
sponse, and neglect correlations between the electrons in 
the two wires. For quantitative comparison with exper- 
iments, one needs to include higher order terms in the 
tunneling and coulomb induced interwire correlations. 
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APPENDIX A: TWO PARTICLES 

The physics of interaction induced electron localization 
is very simply understood by studying two particles in a 
parabolic confinement. The two particle system is partic- 
ularly appealing because it is equivalent to a problem in 
one dimensional single particle quantum mechanics and is 
therefore approachable through elementary means. The 
two and three dimensional version of this problem has 
been extensively studied 18,19 . 

In section IA II we show this exact solution. In sec- 
tion we solve this problem through the Hartree-Fock 
method used in section By comparing these two re- 
sults one quantifies the accuracy of our approximations. 



1. Exact Solution 

In coordinates scaled by the oscillator length w, the 
two particle Schrodinger equation is 



d 2 + d 2 



x\ 



f d {xi - x 2 ) 



ip{x 1 ,x 2 ){Al) 



= Eip(xi,x 2 ), 

where the energy E is scaled by the oscillator energy 
Huj = h 2 /m*w 2 . This equation separates into center of 
mass an relative coordinates A = (x\+x 2 )/2, x — x\ — x 2 
by writing ip(xi,x 2 ) = ^(X)tp(x), which obey 

[(-d 2 x /4 + X 2 ] *(A) = E x <f(X) (A2) 
[-dl + x 2 /A + f d (x)/a] 1>(x) = K1>(x), (A3) 
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FIG. 7: (Color Online) Energy of two electrons in a one di- 
mensional harmonic harmonic well. Thick line: exact singlet; 
thin line: exact triplet; red dotted line: variational singlet; 
blue dashed line: exact Hartree-Fock. Variational Hartree- 
Fock is not shown as it is indistinguishable from the blue 
dashed line. Note logarithmic scale. 




FIG. 8: Relative wavefunction for two particles in a harmonic 
well for (a) a = 1000; (b) a = 10; (c) a = 0.1. All lengths 

are measured in units of d — y^h/mco. Solid line shows exact 
results, while dashed lines show variational approximation. 



with E — E x + Ex- Up to normalization, the lowest 
energy center of mass wavefunction is ^x(X) = e~ x , 
and Ex = 1/2. Our task therefore reduces to solving the 
single particle Schrodinger equation (|A3|I . The electron 
spin enters in that for a spin singlet state ip(x) must be 
even, while for a spin triplet, ip{x) must be odd. 

Figure [7| shows the energy E of the singlet and triplet 
states as a function of a, while figure [H] shows ip(x) for 
several representative values of a. 



take ip{x) oc e~( x ~rf/ 2s2 ± e ~( x ~ t ) 2 / 2s2 where t and s are 
free parameters. The expectation value of the energy is 
then 



g± 




20? 



= f d (x)e- x2/s2 (e tx/s2 ± e~ te/s2 ) (A4) 



For d -C s we extract the leading d dependence of g± to 
arrive at 



dxe x ' 

x\pa 



(e TI - e" TX ) + 0(d 2 /s 2 ) (A5) 



9+ = 9- + 9 

2 r 

9 = -7E h 



(A6) 

7 - log (d 2 /vs 2 )] + 0(d 2 /s 2 ), (A7) 



where r = t/s, 7 « 0.577 is Euler's gamma and v is a 
constant which depends on the exact form of the regu- 
larization. For regularization f^ x \ v = e 2 and for regu- 
larization f^ 2 \ v = 4. The energy is only logarithmically 
sensitive to this constant. 

The separation between the two peaks in the wavefunc- 
tion is given by t, and the width of each peak is given 
by s. These play the role that n _1 and t played in sec- 
tion llVI The lowest energy variational wavefunctions are 
compared to the exact results in figure |S] By expand- 
ing for small t, one sees that for any a < 00, the lowest 
energy variational state always has t =/= 0. 

When appropriately biased, the tunneling current is 
pr oportio nal to \f(k)\ 2 , where f(k) = (l|i/j(fc)|2), k = 
y / 2m*£f—eBd t /Ti, and |1), |2) represent the ground state 
with 1 and 2 particles. Explicitly, 

2 3/4 



/(*) = 



31/2 



dx ip{x)e 



-2ikx/3-x 2 /12-k 2 /6 



(A8) 



Figure El shows this quantity for two illustrative values 
of a. As with the Hartree-Fock theory in section IV Al 
one sees fringes from the fact that each electron is in a 
quantum superposition of two places. The period (in k- 
space) of these fringes becomes smaller as interactions 
increase. 

If one expands the potential in l|A3() about its mini- 
mum, one sees that the size of each peak in the wave- 
function becomes independent of a for sufficiently small 
a. This is born out by analyzing the variational energy 
(|A4|> in the limit of large t. Consequently the envelope 
of the finges in f{k) becomes independent of a for strong 
interactions. 

In a flatter potential (V ~ x 13 with j3 > 2, the peaks 
in ip{x) will become narrower as a increases, resulting in 
a broadening of f(k). 



a. Variational Approximation 



Hartree Fock 



Qualitative features of the exact solution are most 
readily understood through a variational approach. We 



We compare our exact results with the Hartree-Fock 
approximation. For two particles in the same/different 
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FIG. 9: (Color Online) Tunnelling matrix element |/(fc)| 2 for 
two harmonically trapped particles in a singlet state: exact 
a = 100 (solid), a — 0.05 (dashed); hartree fock a = 100 (red 
dotted), a = 0.05 (green dotted, scaled by 0.32). All lengths 
are measured in units of d = ^Jh/mu. Inset shows triplet 
states. 



spin states, Hartree-Fock is equivalent to the variational 
ansatz 

i/j^(x 1 ,x 2 ) = (<Mxi)0 2 (x 2 ) - h(xi)(/)i(x2))Ci^C2^ (A9) 

^n( x i' x 2) = [-M^O-M^Xi^i - 02(xi)0i(x 2 )C2tCh] > 

where <f>i(x) represent single particle orbitals and cor- 
responds to having the j'th particle in the spin up state. 
This ansatz clearly leads to a spin triplet state when the 
two particles have the same spin, but does not lead to any 
particular symmetry when the particles have dissimilar 
spins. 

We focus on the case of dissimilar spins. As might 
be expected from symmetry, we numerically find that 
the minimum energy configuration always has 4>q,{x) = 
4>i{—x) = 4>{x) [which when 4>\(x) is even is equivalent 
to foix) = <t>i(x)]. We then find that the minimal energy 
Hartree-Fock state obeys the self-consistent equation 

d 2 x <j)(x) = {x 2 - nMx) + ?0(ar) [dx \<b(x)\ 2 f d (x + x), 



which is readily solved numerically. 



(A10) 



a. Variational 

As with the exact state, we gain extra insight by a vari- 
ational solution, taking 4>{x) cx e^"*/ 2 ) / s , for which 



(E) 



4 



8 

T 



-t 2 /s 2 



2sa 



9-1 



(AH) 



Not surprisingly, the interaction term in the energy is a 
linear combination of the interaction we previously found 
for triplet and singlet states. 



Unlike when we applied the variational approximation 
to the exact equations (|A4() . here we find that for large a, 
the separation t is strictly zero. As a decreases there is 
a transition where t discontinuously jumps to a nonzero 
value, corresponding to the formation of a spin density 
wave. This discontinuity is also seen in the full numerical 
solution to IjAlOj) . and for d — 0.01 occurs near a = 4.7. 

The variational approximation to the Hartree-Fock en- 
ergy (| Al If) is almost indistinguishable from the exact 
Hartree-Fock energy (|A10() . Figure [7] compares this ap- 
proximate result to the exact energy from solving (|A3|) 
and its variational approximation l|A4|l . As we antici- 
pated, the Hartree-Fock result is nearly exact in the lim- 
its of very strong and very weak interaction. For interme- 
diate coupling however the Hartree-Fock energy deviates 
from the exact result. 

Within Hartree-Fock, the tunneling current is propor- 
tional to \4>k\ 2 , where 4>k — J dxe~ lkx <p(x). Figure 
compares \4>k\ 2 to the exact result for \f(k)\ 2 . Since 
the Hartree-Fock result does not capture the orthogo- 
nality catastrophe physics, the two predictions differ by 
a significant scale factor [0.32 at a = 0.05]. Further- 
more, the exact result contains fringes which are not seen 
in Hartree-Fock. These fringes occur because the exact 
wavefunction places each electron in a superposition of 
two positions. Once rescaled, the Hartree-Fock predic- 
tion gives the envelope of the exact result . 

We expect these observations to be generic. Even in 
the system with more particles, the Hartree-Fock predic- 
tion should match the scaled envelope of the exact result. 
If finite temperature/lifetime effects wash out the fringes 
then the Hartree Fock should correctly predict the depen- 
dance of the tunelling current on momentum (magnetic 
field). 



APPENDIX B: MULTIPLE- WELL POTENTIALS 

The power law potentials used in the body of this paper 
capture the slow variation of electrostatic fields felt by 
electrons in typical experiments. Some features of these 
power-law potentials are overly-simplistic. Most impor- 
tantly, due to experimental constraints on the place- 
ment of gate electrodes, a multiple- well structure may be 
formed along the wire, as illustrated in figure^] Qual- 
itative understanding of the electronic state is found by 
noting that the local chemical potential is small in the 
central region but large near the edges, resulting in a 
low density crystaline state surrounded by a high density 
liquid state. A Hartree-Fock calculation, illustrated in 
figure ITUI verifies this intuitive picture. 

The localized states which are formed above the barrier 
are a purely many-body effect. Similar localized state are 
believed to occur above quantum point contacts^ and 
may give rise to the so called 0.7 anomaljiSi. 
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FIG. 10: (Color Online) Top: Double well potential along 
wire generated by gates. Bottom: Hartree-Fock electron den- 
sity in this potential [see fig. [£] for key] . As illustrated, the 
characteristic length scale of the potential is w. The interac- 
tions are characterized by a = O.lw and d = 0.07«;, yielding 
a low density crystal-like central region surrounded by a high 
density liquid-like region. 
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